Magnetic excitations in multiferroic TbMnOa 



o 
o 

(N 

> 

o 
O 

m 

(N 



i 

-4-* 

S 

i 

O 

o 



> 
o 

(N 
O 

O 

c3 



c 

o 
o 



X 



D. Senff, 1 P. Link, 2 'BK. Hradil, 3 A. Hiess, 4 L.P. Regnault, 5 Y. Sidis, 6 N. Aliouane, 7 D.N. Argyriou, 7 and M. Braden 1 '! 

II. Physikalisches Institut, Universitat zu Koln, Ziilpicher Str. 77, D-50937 Koln, Germany 
2 Forschungsneutronenquelle Heinz Maier- Leibnitz (FRM II), 
TU Miinchen, Lichtenbergstr. 1, 85747 Garching , Germany 
3 Institut fur Physikalische Chemie, Universitat Gottingen, Tammanstr. 6, 37077 Gottingen , Germany 
^Institut Laue-Langevin, Boite Postale 156, 38042 Grenoble Cedex 9, France 
5 CEA-Grenoble, DRFMC-SPSMS-MDN, 17 rue des Martyrs, 38054 Grenoble Cedex 9, France 
6 Laboratoire Leon Brillouin, C.E.A./C.N.R.S., F-91191 Gif-sur-Yvette Cedex, France 
7 Hahn-Meitner-Institut, Glienicker Str. 100, 14109 Berlin, Germany 
(Dated: February 6, 2008, preprint) 

The magnetic excitations in multiferroic TbMnOs have been studied by inelastic neutron scatter- 
ing in the spiral and sinusoidally ordered phases. At the incommensurate magnetic zone center of 
the spiral phase, we find three low-lying magnons whose character has been fully determined using 
neutron-polarization analysis. The excitation at the lowest energy is the sliding mode of the spiral, 
and two modes at 1.1 and 2.5 meV correspond to rotations of the spiral rotation plane. These latter 
modes are expected to couple to the electric polarization. The 2.5 meV-mode is in perfect agreement 
with recent infra-red-spectroscopy data giving strong support to its interpretation as an hybridized 
phonon-magnon excitation. 

PACS numbers: 75.30.Ds, 75.47.Lx, 75.40.Gb, 75.50.Ee, 77.80.-e 



The perovskite TbMn0 3 [l| is a key material for the 
new class of multiferroic transition metal oxides 0, 0, , 
since the polarization in the ferroelectric phase is size- 
able and since the magnetoelectric coupling is remark- 
ably large 0, 0. The antiferro-type ordering of the 
single-occupied e g -orbitals implies a ferromagnetic near- 
est neighbor (nn) interaction J FM in LaMn0 3 as well 
as in TbMn0 3 [6|, but in TbMn0 3 , a large c-axis oc- 
tahedron rotation [7j yields a sizeable overlap of the eg- 
orbitals of next-nearest neighbor sites (nnn) along b ren- 
dering the associated magnetic interaction strongly anti- 
ferromagnetic, J nnn The well-defined ferromagnetic 
order in the a, 6 planes of LaMn0 3 , therefore, becomes 
strongly frustrated in TbMn0 3 giving rise to the complex 
magnetic ordering which finally causes the multiferroic 
behavior. 

At T N =42 K, the Mn spins in TbMn0 3 order in a lon- 
gitudinal spin-density wave (SDW) with a wave-vector 
of q=(0, 0.28,0) in reduced lattice units of the Pbnm- 
structure Upon further cooling, the modulation 

vector changes slightly until at Tn=28 K a second transi- 
tion into a spiral phase occurs [jj, |£| • Here, the magnetic 
order corresponds to an elliptic cycloid still modulated 
along the fe-direction but with the spins rotating around 
the a-axis . Associated with this SDW to spiral transi- 
tion, the spontaneous electric polarization parallel to the 
c-axis appears. 

The appearance of the electric polarization at the tran- 
sition into the spiral phase was microscopically explained 
by Katsura et al. PJJ]. This idea was followed by studies 
of the phenomenology an d its extension to the lattice 
interaction Due to the non-collinear spin arrange- 

ment in the spiral phase, the inverse of the Dzyaloshinski- 
Moriya interaction implies an uniform displacement with 
an electric polarization given by: 
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FIG. 1: (color online) Raw-data scans at the incommensurate 
zone center and at the 6-axis zone boundary in the spiral 
phase (PANDA and 4F, fc/^l^A" 1 ) a); Q-dependence of the 
signal b), the insert c) compares the suppression of the 4.5 
meV feature with the expectation due to a Tb magnetic form 
factor (PUMA, fci=2.66A~ 1 ). Full polarization analysis of the 
magnetic excitations at the incommensurate zone center and 
at the zone boundary (IN14, = 5A _1 ) d) and e). 



Pocry x (S; x Sj) (1), 

for a pair of spins Si,Sj with a distance vector ry. Con- 
sidering the mechanism in equation (1) there is a close 
coupling between the dielectric properties and the mag- 
netic excitations which should lead to hybridized phonon- 
magnon excitations. First evidence for such mixed exci- 
tations was recently obtained in Infra-Red (IR) optical- 
spectroscopy measurements [l3| .however, their role in the 
multiferroic behavior is unclear as their magnetic coun- 
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FIG. 2: Raw-data scans to determine the magnon dispersion 
in TbMn0 3 along a and c directions starting at the incom- 
mensurate zone center (0,0.28,1), (PUMA, fc l =2.662A _1 ). 

terparts have been little explored thus far [T^ . 

In this letter we present the results of inelastic neutron 
scattering (INS) experiments to unravel the magnetic ex- 
citations in TbMn0 3 , in particular the low-energy modes 
relevant for the multiferroic behavior whose characters 
were determined using longitudinal polarization analysis. 
The excellent agreement with the IR studies strengthens 
the evidence for mixed excitations. 

A ~900 mm 3 large untwinned single crystal of 
TbMn0 3 was grown using an image furnace hj. The 
crystal exhibits the magnetic transitions at 42 and 28 K 
in agreement with previous studies Q, Q . INS experi- 
ments were performed on various triple-axis spectrome- 
ters : on 4F (cold neutrons, Laboratoire Leon Brillouin), 
on PANDA (cold neutrons) and PUMA (thermal beam, 
both at the Forschungsreaktor Munchen II) and using full 
polarization analysis with the CRYOPAD setup on IN14 
(cold neutrons, Institut Laue Langevin). Measurements 
were taken in the ferroelectric (T=17 and 23 K), in the 
paraelectric SDW and in the paramagnetic phases. 

Fig. 1 shows scans taken at the incommensurate zone 
center Q=(0, 0.28,1) and at the zone boundary in b di- 
rection in the spiral phase at T=17 K. Typical scans to 
determine the dispersion along the a and c-directions are 
shown in Fig. 2; along these directions the frequencies 
rapidly increase. In all scans we find an excitation at 
4.5 meV which in agreement with the interpretation in 
reference 0] can be identified as a Tb crystal-field exci- 
tation, since it does not disperse and since it exhibits a 
weaker dependence on the length of the scattering vec- 
tor, |Q|, compared to the other contributions. The | (In- 
dependence, furthermore, is in good agreement with the 
Tb magnetic form factor, see Fig. 1 c). 

Fig. 3 shows the magnon dispersion of TbMn0 3 
obtained in the spiral (T=23 and 17 K) and in the 
paraelectric SDW phase (T=32 K). The dotted lines 
in Fig. 3 denote the magnon dispersion measured 
in LaMn0 3 0|. Only along the antiferromagnetic 
c-direction, LaMn0 3 and TbMn0 3 exhibit a qualita- 
tively similar dispersion. Along the a, b plane, the dis- 
persion is essentially flattened in TbMn0 3 . The up- 
per cut-off energy in TbMn0 3 amounts to ^8 meV, 
whereas magnon energies in LaMn0 3 extend up to 33 
meV [lij . The pronounced softening of the magnon fre- 



quencies for wave vectors along the ferromagnetic planes 
reflects the frustrating nnn interactions and the weak- 
ened nearest neighbor interaction 0. The magnons in 
TbMn0 3 are broadened, and we find evidence for a split- 
ting of excitations even in the SDW phase, which, how- 
ever, requires further studies. Tentatively, we may de- 
scribe the dispersion along a and c in a S=2 Heisenberg 
model with a single magnon branch including single-ion 
anisotropy, D. Thereby, we obtain the interaction param- 
eters J FM =0.15(l)meV (nearest neighbors along the a, b 
planes) and J AFM =-0.31(2)meV (neighbors along c) and 
D=0.09(l)meV. These values should be compared to the 
values of J PM =0.83 meV and Jafm~ -0-58 meV observed 
in LaMn0 3 [lfj . In the spiral phase, mode splitting 
is observed throughout the Brillouin zone and a more 
sophisticated model to fully describe the spin-excitation 
dispersion is needed. 

Concerning the high-energy part of the dispersion, our 
data agree qualitatively with the study of the magnon 
dispersion by Kajimoto et al. Q, but the splitting in the 
excitations has not been analyzed in reference [14| prob- 
ably due to the fact that this group did not analyze the 
dispersion starting at the incommensurate zone center. 
In addition, we have explored in detail the low-energy 
part of the dispersion at the incommensurate zone cen- 
ter and the soft branches along the 6-direction. As seen in 
Fig. 1 and in the dispersion along the b direction shown 
in Fig. 3 a) , there are actually three low-energy branches 
emerging out of Q=(0, 0.28,1) . In view of the magneto- 
electric effect, these are the most relevant modes. The 
unpolarized data in Fig. 1 a) show the three contribu- 
tions (plus the Tb crystal field at 4.5 meV), but this type 
of data does not allow us to determine unambiguously the 
polarization of the different contributions in spite of the 
many different Brillouin zones explored. 

We have employed longitudinal polarization analysis 
yielding important additional information : in the spin- 
flip scattering channel only those excitations are de- 
tected whose polarization is perpendicular to the chosen 
neutron-polarization axis |17j . We have performed the 
experiment with the a direction vertical to the scattering 
plane using the CRYOPAD-device on IN14. The setup 
yielded a remarkable precision with a flipping ratio of 

JNSF . jSF = 35 and 

an accuracy of about 2% in the 
transverse polarization terms. With the three neutron- 
polarization directions : parallel to Q, (x), parallel to a, 
(z), and perpendicular to both Q and a, (y), one may 
distinguish between excitations polarized perpendicular 
to the spiral plane in TbMn0 3 (i.e. perpendicular to the 
b, c spiral plane, S±) and those polarized within the spi- 
ral plane (i.e. in the b, c plane, Su). In the spin-flip (x)- 
channel the sum of both magnetic contributions appear, 
whereas the spin-flip (y) and (z) channels measure only 
the S± and the S\\ contribution, respectively ^t|- By 
subtracting the intensities obtained in different channels 
one directly obtains S± and S\\ without any background 
assumption. The two contributions S± and 6*11 are sep- 
arated in the lower part of Fig. 1. At the incommensu- 
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FIG. 3: Dispersion of the spin-wave excitations in TbMn0 3 
along the a and the c directions of the Pbnm lattice in the 
paraelectric sinusoidal phase a) and in the ferroelectric spiral 
phase b) . Dispersion along the modulation vector in the 
spiral and in the SDW phase c). Polarization of the three 
low-energy modes at the incommensurate zone center d) : 
the two a-polarized modes, (i) and (ii), and the sliding mode 
of the spiral, (iii), see text. 



rate zone center, the excitation with the lowest energy 
of 0.2(1) meV is b, c polarized. Since one may exclude a 
longitudinal spin excitation in this high-moment ordered 
structure, this mode must be attributed to the trans- 
verse magnon polarized in the spiral plane. This mode 
is the sliding mode of the modulated magnetic structure 
its polarization scheme is shown in Fig 3d) (iii). It has 
no impact on equation (1) describing the magnetoelec- 
tric coupling. Therefore, it should not be relevant for the 
multiferroic behavior. Since the spiral in TbMn0 3 is not 
perfectly circular an anisotropy term can explain a finite 
energy of the sliding mode in addition to pinning effects. 

In the S± channel, there are two contributions at ener- 
gies of 1.08(8) and 2.46(9) meV. The low-energy excita- 
tion of a multiferroic system where fcrroclcctricity origi- 
nates from the inverse Dzyaloshinski Moriya interaction 
has been discussed recently by Katsura et al. |l8|. Be- 
sides the sliding mode and the phonons, there are indeed 
two magnetic excitations with a spin polarization per- 
pendicular to the spiral plane (o, c plane in the case of 
TbMn0 3 ) corresponding to S±. These modes are illus- 
trated in Fig. 3d) i) and ii). For simplicity, we assume 
that the phase shift between the b and the c components 
of the spiral in TbMn0 3 is exactly 7r/2 and that the two 
coefficients are of equal size, i.e. the spiral is circular. 
Then, one may describe the magnetic order by : 



Si = S ■ cos(q spil . a i • Ri) • e y + S • sin(q spira i • R ; ) • e z (2), 



where Ri is the position vector of the spin Si. An 
a-polarized magnetic excitation at the incommensurate 
zone center necessarily possesses the same modulation 
along 6, however the oscillating part can be in-phase ei- 
ther with the b component (cosine term in equation (2), 
see Fig. 3d) (ii) ) or with the c component (sine term in 
equation (2), see Fig. 3d) (i) ). The arising polarization 
patterns are distinct : in the first case the spiral plane 
rotates around the c direction and it rotates around b 
in the second case. Correspondingly, the cross product 
Si x Sj between neighboring spins rotates around the c 
direction and around the b direction in these cases. Since 
the modulation vector or ry in equation (1) remain al- 
ways parallel b, the rotation of the Si x Sj cross product 
through the magnon has a strong linear coupling to the 
polarization in the second case, which we label w_ fol- 
lowing reference ^3- The w_-excitation is the Goldstone 
boson of the multiferroic phase, it may shift to finite en- 
ergy due to anisotropy terms. It is a magnon-phonon 
hybridization with the phonon describing the electronic 
polarization along the a direction. The other transverse 
magnon couples only quadratically, and thus weakly, to 
the polarization. It weakly modulates the static polariza- 
tion parallel to c. Both magnons posses finite frequencies 
due to anisotropy effects. 

Recently, Pimenov et al. EH have reported the ob- 
servation of a low-energy excitation in IR spectroscopy 
studies of GdMn0 3 and TbMn0 3 which could be fully 
suppressed by applying an external magnetic field. They 
interpreted this excitation as an electromagnon, the hy- 
bridized magnetic and phononic excitation. This result 
together with the analysis given in reference |l8| agree 
perfectly with our INS study of the magnetic excitations. 
The electromagnon feature seen by Pimenov et al. at 
20cm _1 =2.48 meV agrees with one of the low-lying mag- 
netic excitations in the S± contribution. We, therefore, 
tentatively attribute the 2.46 meV neutron-scattering 
excitation to the lj_ Goldstone boson. According to the 
analysis given above, this mode should couple to the po- 
larization in the a-direction, which is indeed the direction 
of electric field in the IR experiment sensing the lattice 
part of the excitation. In TbMn0 3 the spiral order might 
be more complex than the ideal circular cycloid consid- 
ered in equation (2); therefore the characters of the two 
S± -polarized magnetic excitations might weakly mix ren- 
dering both IR-active. It is interesting to note that also 
the energy of the lower S± mode agrees perfectly with 
the - though weaker - feature in the IR response at 1.24 
meV in reference E3l- We emphasize, that associating 
the 2.46 meV neutron-scattering excitation with the uj~- 
mode only bases on the comparison with the IR results; 
actually one might expect this mode to exhibit a lower 
energy since it describes the flip of the spiral into the 
a, 6-plane which seems to occur in TbMn0 3 under field 
as well as in related ReMn0 3 US El. 

Along the b direction the three low-energy branches 
can be followed till the incommensurate zone bound- 
ary and even beyond in the spiral phase. The disper- 
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FIG. 4: (color online) Temperature dependence of the mag- 
netic scattering at the incommensurate zone center deter- 
mined with unpolarized neutrons (PUMA, fci=2.662 A -1 con- 
stant) a) and with polarized neutrons (IN14, fc/ = 1.5 A -1 con- 
stant) c); the inset b) shows the magnetic elastic scattering 
at Q=(0, 1.72,1) as function of temperature. 

sion unambiguously shows that anisotropy terms play 
a dominant role in TbMn0 3 ; they appear to determine 
the sequence of magnetic phases as function of temper- 
ature and magnetic field. Introducing the frustrating 
J nnn interaction into the non-frustrated magnon disper- 
sion along b, see insert in Fig. 3a), one may understand 
the 6-dispersion in the incommensurate phase. The frus- 
tration generates a dispersion minimum whose position 
as well as the analysis of the exchange energy in the or- 
dered state allow one to estimate the interaction ratio 
to be of the order of J nnn / J FM =0.8. The finite ener- 
gies of both the feature and the other a-polarized 
zone-center magnetic excitation arise from either a single- 
ion anisotropy, from exchange anisotropy or from the 
electron-lattice coupling. More experimental and theo- 
retical efforts are, however, needed in order to fully char- 
acterize the aspects arising from folding the magnon dis- 
persion into the incommensurate ordering. 

The inset of Fig. 4 shows the spin-flip signal at the 
(0,1.72,1) magnetic Bragg reflection which essentially de- 
pends on the c-component of the ordered moments; its 
temperature dependence clearly shows the drastic in- 
crease of the c-component in the spiral phase j§|. The 
temperature dependence of the spectra obtained at the 
incommensurate zone center, see Fig. 4, show that the 
two well-separated a-polarized excitations in the spiral 
phase merge into a broad signal upon heating into the 
SDW phase. In analogy to the discussion of the low- 
lying excitations in the spiral phase, one may expect sev- 
eral contributions in the SDW phase as well. Single-ion 
anisotropy causes a finite energy for the transverse mag- 
netic excitations and will yield a splitting for polarization 
parallel a and c. Furthermore, there will be a separation 
due to the phase of the oscillating part which may be 
either in-phase with the SDW or shifted by ir/2. The in- 
phase modes are irrelevant for the magnetoelectric cou- 
pling at least within the Dzyaloshinski-Moriya scenario of 
equation (1) since they lead to a collinear spin structure. 
In contrast, the two 7r/2-phase modes may couple with 



electric polarization even in the paraelectric SDW phase 
for both magnetic polarizations. Besides the broaden- 
ing of the neutron signal in the SDW phase, there is a 
temperature driven change in the spectral weights for 
the different polarizations. There is still a low-lying S\\ 
contribution in the SDW phase in agreement with the 
character of the transition into the spiral phase where or- 
dered moments are displaced parallel to the c direction. 
This low-lying S\\ signal can be considered as the quasi- 
soft mode associated with the SDW to spiral transition. 
This mode should be strongly coupled to the polariza- 
tion along c. The total spectral weight of the a- or 
polarized contributions, however, seems not to be focused 
in a single mode but extends to higher energies. Compar- 
ing the response in the Sj_-channel in both phases, one 
may conclude that the w_ mode is little changing across 
the transition, whereas the other S^-contribution softens 
in the SDW phase. These findings agree once more with 
the IR studies, where the electromagnon for fields along 
a persists into the paraelectric SDW phase. 

In conclusion the INS studies of the magnetic exci- 
tations in multifcrroic TbMn0 3 with and without polar- 
ization analysis allow us to identify the frequencies and 
character of the low-energy modes associated with the 
magnetoelectric coupling. The good agreement of the 
observed frequencies with those of recent IR studies give 
strong support to the interpretation as a magnon-phonon 
hybridized electromagnon mode. 
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